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\^ • Gamow-Teller (GT) resonances in finite nuclei are studied in a fully consistent 

relativistic random phase approximation (RPA) framework. A relativistic form of 
the Landau-Migdal contact interaction in the spin-isospin channel is adopted. This 
' choice ensures that the GT excitation energy in nuclear matter is correctly repro- 

00 ; 

■ duced in the non-relativistic limit. The GT response functions of doubly magic nuclei 

^^Ca, ^''Zr and ^osp]^ calculated using the parameter set NL3 and (7g=0.6 . It is 



X 



found that effects related to Dirac sea states account for a reduction of 6-7% in the 
GT sum rule. 

PACS numbers: 21.60.-n, 21.60.Jz, 24.10.Jv, 24.30. Cz 

In recent years, the relativistic mean field (RMF) theory with non-linear meson self- 
interactions has achieved a great success in describing bulk properties of nuclei, not only 
spherical but also deformed nuclei and nuclei far from the /5-stability lineQ]- In particular, 
a fully consistent relativistic random phase approximation (RPA) based on the RMF has been 
established . The consistency implies that the particle-hole (p-h) residual interaction 

and the nuclear mean field are calculated from the same effective Lagrangian. The relativistic 
RPA is equivalent to the ti..depe„de„t RMF the s,„all a,„p.itude l.^itfl Q only ,f the 
particle-hole configuration space includes not only the pairs formed from the occupied and unoc- 
cupied Fermi states but also the pairs formed from the Dirac states and occupied Fermi states^]. 
It has been found that the effective Lagrangians which can well describe the ground state prop- 
erties of nuclei could also reproduce their collective excited states and giant resonances^]. 

The investigation of nuclear excitations involving spin-isospin Gamow-Teller(GT) resonances 
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has attracted a great interests for a long time. The nuclear GT transitions are essential to many 
important processes in particle physics and astrophysics related to neutrino-nucleus interactions. 
The quenching of the observed GT strength in nuclei is a long standing issue both experimen- 
tally and theoretically. It has been advanced through the recent extensive investigations of GT 



resonances by the charge exchange reactions at intermediate energies |10|Jll|. The recent {p,n) 
and (n, p) measurements have observed about 90% of the Ikeda sum rule|l2l| at energies < 50 
MeV in medium heavy nuclei Theoretical investigations have shown that there are several 
possible mechanisms for the unseen strength, which range from the admixture of 2p — 2h compo- 
nents in even-even nuclei^^ to the coupling to the A— hole sector^]. The theoretical study of 
GT matrix elements in nuclei has already revealed some valuable information of nuclear struc- 
ture effects and non-nucleonic degrees of freedom such as mesons and delta-isobars. A further 
and clear understanding of the GT resonances is still required. Up to now most investigations in 
the relativistic approach are restricted to the electric excitations. Although the GT resonance 
has been recently studied in the relativistic RPA 3, Q| a proper relativistic form of the p-h 
interaction in the spin-isospin channel has not been obtained. 

A relativistic form of the Landau-Migdal g' term which could reproduce the non-relativistic 
result for the GT excitation energy in nuclear matter has been recently proposed 17] . The study 
of the relativistic description of GT resonances in nuclear matter has pointed out a new quenching 
mechanism due to the relativistic effects. The purpose of this letter is to investigate the GT 
resonances in finite nuclei in the fully consistent relativistic RPA. The GT response functions of 
^°^Pb, ^°Zr and ^'^Ca are studied. We shall emphasize the contribution of the Dirac states to the 
quenching of the Ikeda sum rule in the GT stren gth s due to the completeness of nuclear wave 
functions as pointed out for nuclear matter in ref. 17 1. 

The relativistic RPA is built on the ground state which is calculated in the RMF theory. For 
a consistent RPA calculation the residual p-h interaction giving rise to giant resonances must be 
obtained from the same Lagrangian. The spin-isospin correlations in the relativistic approach 
are induced by the isovector mesons vr and p. In the existing effective Lagrangians explicitly 
built for RMF the pion term does not appear since it does not contribute to the mean field 
while for the p the tensor coupling term is generally ignored and the vector coupling term is 
strongly renormalized to give the correct neutron-proton symmetry energy. To remain consistent 
we choose to introduce in the effective Lagrangian a pseudo-vector pion term but no p tensor 
term. Thus, the GT correlations will be induced by the Lagrangian 

fn— — 

C = -ip-f^l^d^T ■ n^fj - gptp-ff^df^T ■ p^ip . (1) 
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It is well known that due to the pseudo- vector coupling of the pion and tensor coupling of the 
p an effective interaction with the Landau-Migdal parameter g' term has to be introduced, which 
has the form g'cri ■ (T2T1 ■ T2 in the non-relativistic limit. There are several ways to introduce 
this term in the interaction LagrangianQ, Q]. However, It has been pointed out in R,ef . ^\ 
that the way to introduce g' in the relativistic model is model-dependent and that the following 
choice 

£ = ^^rf ^V^r^^^A, rf = 757^r„ g, = g'^, (2) 

insures that, in the non-relativistic limit one can recover the expression of the GT excitation 
energy in nuclear matter. 

The response function of the system to an external field is given by the imaginary part of the 
retarded polarization operator, 

i?(P, P; k'; E) = -Imn^(P, P; k, k'; E) , (3) 

TT 

where P is the external field operator. The relativistic RPA polarization operator is obtained 
by solving the Bethe-Salpeter equation jiol. 

n(Q, Q'- k, k', E) = no(Q, Q'; k, k', E) 
[ d'k,d%Uo{Q,r;k,k,,E) 

A(ki,k2,E)n(ri,Q';k2,k',P) , (4) 

where the Fj's and Z^j's are the vertex couplings and the corresponding meson propagators, 
respectively, Q and Q' represent the external or vertex operators. Generally, the operators of 
spin-isospin excitations are 

1 ^ 

* i 

where t± = (r^ ± iry) / and L=0 for the GT operator. The operator is usually multiplied by 
7o for vector density excitations. 

The mo del- independent Ikeda sum rulej3] is expressed as 

(0|P+P_|0) - (0|P_P+|0) = 3(iV - Z) , (6) 

where |0) is either the nuclear correlated or uncorrelated ground state. 
The vertex operators corresponding to vr- and p-coupling are 

r.r^ = 757Mg'^(l)757/.g^(2)riT2 , 
T,,T>; = 7m(1)7'^TiT2 

= [70(1)70(2) - 7o750-(l)7o75cr(2)]TiT2 , (7) 
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with the corresponding meson propagators 
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{27r)^ — + irj 
1 



(27r)'^ — m? + irj 



(9) 



The vertex operator of the Landau-Migdal force introduced in Eq.Q is 



r^r'' = [7075(1)7075(2) - 7ocr(l)7ocr(2)]TiT2 



(10) 



The time component and the current part have to be calculated separately. The propagator for 
the contact interaction term in the Lagrangian is 

2 
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The multipole expansion in momentum space is performed and the integrals of the angular 
parts can be carried out. The retarded unperturbed polarization operator with a fixed angular 
momentum J in the momentum space is written as 



U',{Q,Q';k,q;E) 




{Sh - Bp) + E + irj 



{Sh - Ep) - E - irj 



ah 



{sa - Eh) + E + ir] 



(12) 



[Ea - Eh) - E - ir] 

where h, p and d correspond to the occupied states in the Fermi sea, positive energy unoccupied 
states and negative energy states in the Dirac sea, respectively. The advantage of solving the 
Bethe-Salpeter equation in the momentum space is to allow the inclusion of as many configu- 
rations as required. Indeed, the numerical work amounts to invert matrices the size of which 
depends on the number of mesh points chosen in momentum space and the number or vertex 
operators. Therefore, increasing the number of configurations does not affect substantially the 
numerical effort in contrast to the matrix diagonalization method. 

The external operator V after multipole expansion can be expressed as 



P = 7or^-'[K 



j-i a J . 



(13) 
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The vertex operator produced by the pion is 



L=J±1 

L 1 




lo[YL(^(r]jjL{kr) . (14) 



The time component of the operator produced by the vector part of the p vanishes and only its 
space component contributes and leaves a minus sign on the propagator: 

Q' = lol5[Yj(E)a]jjjikr) . (15) 

The time component of the vertex operator produced by the Landau-Migdal force is: 

Q'o = 7ol5Yj{r)jj{kr) , (16) 

Its space components are: 

= 70 [Yj±i ® cT]jjj±^{kr) , (17) 

which produces also a minus sign on the propagator. The isospin part of the operators is 
multiplied by r_|_/-\/2- There are in total 5 operators for Q, therefore the polarization 11'^ {Q, Q') 
is a 5 X 5 matrix. 

The single-particle energies and wave functions are solutions of the self-consistent RMF equa- 
tions obtained with the starting effective Lagrangian. A nucleon state can be specified by a set 
of quantum numbers a = {uaja, jayf^a) = (o-, "^a), where ga=-l and +1 for neutron and proton 
states, respectively. The nucleon wave function with energy is written as 

V'^ = - I '^"^""^ I yair)x{<la) ■ {li 

Fjr)(T ■ r 



xila) and being the isospinor and spin-spherical harmonic function, respectively. 

The charge-exchange excitations flip the isospin and pick up only neutron-proton or proton- 
neutron pairs. The reduced matrix element for the external spin-isospin excitation is 



{-^hWVUv) = J {GHG,{h\\Tj., j\\p) 

+F^F,{h\\Tj-i j\\v)y-'dr , (19) 

where h = Uhlhjh, and / is the orbital angular momentum of the lower component in the Dirac 
spinor. Here, we Win.od.ced(.||r„||;)H..U.II|.-.«.l.llV*)Q/T.e .educed n,.^^^ 
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elements of the operators ()14II17|) are 



f(L) 

{nhhjh\\Yj\\nplpjp) + q ^ , (20) 

L=j±i V 2^ + ^ 



X J {GhGp{h\\TLj\\p) + Fj,Fp{h\\TLj\\p)jL{kr)dr 
where f{L) = J+ lforL = J+ l and J for L = J - 1, 



{^h\\Q'W = j {GnFp{h\\Tjj\\p) 

-FnGM\Tjj\\v))3jdr , (21) 

ii^hWQoW = j{GhFp-F,Gp)jj{kr)dr 

{nhlhjh\\Yj\\nplpjp) , (22) 



{4^Hm±iW = J iG,Gp{h\\Tj^, j\\p) 

+FhFpCh\\Tj^, j\\p))jj±iikr)dr . (23) 

We calculate the L=0 GT excitations for the double closed shell nuclei ^^Ca, ^^Zt and ^°^Pb, 
which have been extensively investigated experimentally. The ground state wave functions of 
those nuclei are calculated in the RMF with the parameter set NL3|7|. The continuous single- 
particle spectrum is discretized in a harmonic oscillator basis. An averaging parameter A = 2 
MeV is used to smooth out the response functions, which is performed by replacing the excitation 
energy E + if] in Eq.® by E + iA/2. The residual p-h interaction is produced through the 
isovector mesons and the Landau-Migdal force in the Lagrangian Eqs.(PI2I), where g'=0.6 and 
the standard values /^/47r = 0.08, m.„ = 138 MeV are adopted. All other coupling constants are 
those of NL3. 

The calculated Hartree and correlated GT strengths for ^^Ca, ^"Zr and ^°^Pb are shown in 
Fig.l. The dashed and solid curves correspond to the Hartree and RPA results, respectively. 
The main contributions to the GT excitation are from the neutrons at the Fermi surface excited 
to the unoccupied proton states with the same orbital angular momentum. For instance, the 
main p-h configurations are: {7!'lgQ/2{i^lg9/2)~^) with Eph = 6.3 MeV and (vrl(77/2(z^l5'9/2)~^) 
with Sph = 13.9 MeV in ^^Zr; (7rl/7/2(z/l/7/2)~'), Sp^ = 0.3 MeV and (711/5/2(^/1/7/2)-'), Sp^ = 
8.4 MeV in ^^Ca, which exhibit peaks in the Hartree strengths. There are more contributing 
configurations in ^°'^Pb. The collectivity effect is mainly due to the Landau-Migdal term of the 
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E ( MeV ) 



FIG. 1: GT response functions for ^^Ca, ^"Zr and ^"*^Pb. The GT strengths are calculated in the 
relativistic RPA with the parameter set NL3. The solid and dashed curves correspond to the RPA and 
Hartree strengths, respectively. Experimental energies are shown by the arrows. 

residual interaction. It results in an upward shift of the strength. The centroid energies of the 
main peaks in Hartree and in RPA are listed in Table 1. A good agreement with the experimental 
data is observed when g' = 0.6 is chosen. The total strengths up to 60 MeV exhaust the Ikeda 
sum rule by about 92-94 %. From Table 1 it can also be seen that the total GT strength below 
60 MeV is similar in Hartree and in RPA. In ref.jlTl it was also found that the GT quenching 
is about 12% in infinite matter and 6% in nuclei calculated in Hartree approximation. The 
rest of the strength is connected with the effects of the Dirac sea states. NN pairs involving 
states in the Fermi sea and negative energy states can also carry GT strength, if one considers 
backward-going graphs and the no-sea approximation. 

In summary, we have investigated in finite nuclei the nuclear spin-isospin excitations, espe- 
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TABLE I: GT energies of ^^Ca, ^^Zr and ^^Spb. E corresponds to the centroid energies of theoretical 
results with g' = 0.6. The percentages of Ikeda sum rule are calculated up to < 60 MeV. The 



experimental data are taken from Refs. 



Id, 



21, 



Theoretical Results Experiment 



Nuclei E (MeV) Strength 


% of sum rule 


Energy 


Strength 


RRPA 


Hartree 


RRPA 


(MeV) 


% of sum rule 


48Ca 10.1 


93.0 


93.8 


10.5 


35 


90Zr 15.4 


92.3 


93.2 


15.6±0.3 


90±5 


208pb 18.9 


91.6 


92.6 


19.2±0.2 


60~70 



daily the GT resonance in a fully relativistic RPA. A relativistic form of the Landau-Migdal 
contact interaction in the spin-isospin channel is adopted, which reproduces in the non-relativistic 



limit the excitation energy of the giant GT resonance in infinite matter [17[. The GT resonances 
in ^^Ca, ^°Zr and ^osp]^ investigated. It is found that the RPA strengths up to 60 MeV 
exhaust the Ikeda sum rule by only about 93%. The missing fraction is taken by the pairs 
formed between states in the Fermi sea and Dirac sea. This quenching mechanism is specific of 
the relativistic description of the GT mode. 
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